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The dynamics of dislocations confined to grain boundaries in a striped system are studied using 
electroconvection in the nematic liquid crystal N4. In electroconvection, a striped pattern of con- 
vection rolls forms for sufficiently high driving voltages. We consider the case of a rapid change in 
the voltage that takes the system from a uniform state to a state consisting of striped domains with 
two different wavevectors. The domains are separated by domain walls along one axis and a grain 
boundary of dislocations in the perpendicular direction. The pattern evolves through dislocation 
motion parallel to the domain walls. We report on features of the dislocation dynamics. The kinetics 
of the domain motion are quantified using three measures: dislocation density, average domain wall 
length, and the total domain wall length per area. All three quantities exhibit behavior consistent 
with power law evolution in time, with the defect density decaying as t~^^^, the average domain 
wall length growing as t^^^, and the total domain wall length decaying as t~^^^. The two different 
exponents are indicative of the anisotropic growth of domains in the system. 

PACS numbers: 47.54.+r, 64.60.Cn 



I. INTRODUCTION 

The dynamics of topological defects are observed to 
dominate the temporal evolution of patterns in many 
physical systems. However, our understanding of the 
quantitative contribution of the defect dynamics to the 
evolution of patterns is still not complete. One area in 
which topological defects potentially play a central role 
is the growth of domains in a patterned system after a 
sudden change in the external parameters. For many 
striped systems, such as diblock copolymers and convec- 
tion in fluids, many of the topological defects (such as 
disclinations and domain walls) are the same as those 
that exist in classic models for the growth of "uniform 
domains" , such as Ising models or x-y models For uni- 
form systems, the contribution of the topological defects 
to the time evolution of the domains is relatively well 
understood 0; whereas, this is not the case for striped 
systems. 

It is useful to briefly review the situation for uniform 
domain growth 0. One is generally interested in the 
evolution of a system after a rapid change of an exter- 
nal parameter, or a quench. Typically, one considers an 
initially uniform state that immediately after the quench 
is no longer an equilibrium or steady-state phase of the 
system. The new equilibrium phase is degenerate, and 
domains of the different states form. For example, in an 
Ising system, one would have domains of up and down 
spins. During the subsequent evolution of the system, 
or coarsening, the domains are characterized by a single 
length scale. This length grows as a power of time, i", 
where n is the growth exponent. We use the designation 
"uniform domains" to refer to systems in which within 
a domain, the system is uniform. Examples of this type 
of system include magnetic systems, metallic alloys, bi- 
nary fluids, and nematic liquid crystals One of the 
main goals of this field is to understand the possible val- 
ues of n and other features of the late time scaling. For 



uniform systems, much of the late-time behavior can be 
understood in terms of the topological defects of the or- 
der parameter 0- For most cases, the growth exponent 
is 1/3 if the order parameter is conserved, such as in bi- 
nary fluids, and the exponent is 1/2 if it is not conserved, 
as in magnetic systems 0. Various interesting alternate 
cases exist, especially for the x-y model in one and two 
dimensions 0- In contrast, no similar general framework 
exists for patterned domains. These are systems in which 
within a domain, the system exhibits a pattern, such as 
stripes. 

Stripes, or more generally patterns, occur in a wide 
range of systems including convecting fluids, an- 

imal coats, polymer melts, and ferromagnets. Stripes 
occur both as an equilibrium state of the system, such as 
in diblock copolymers, and as a result of external driving, 
as in convection in fluids. A sudden change of an exter- 
nal parameter, a quench, can bring the system from a 
spatially uniform state to a striped state that undergoes 
coarsening, or phase ordering. The coarsening of striped 
domains has focused on isotropic systems 0, |a IS S 11101 j 
with a significant fraction of the work focusing on simu- 
lations of model equations, such as the Swift-Hohenberg 
equation 0, |H i£i M)] ■ Experimental work has been 
limited to two systems, one isotropic (diblock copolymers 
0,1^) and one anisotropic (electroconvection in nematic 
liquid crystals J2J). For isotropic systems, the stripes 
can have any orientation. In this case, the dominant de- 
fects are disclinations and domain walls, though disloca- 
tions are observed as well. For diblock copolymers, the 
measured growth exponent is consistent with n = 1/4. 
Detailed studies of the disclination dynamics are able to 
explain this measured exponent 

For simulations of isotropic systems, growth exponents 
are usually consistent with 1/4 l3,'5','^ or 1/5 0,1111 [ill. 
Factors affecting the measured exponent include external 
noise and the quantity used to characterize the domains. 
There is evidence that for small enough quench depths 
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the growth exponent is 1/3 An open ques- 

tion is the connection between the simulations and the 
experiments with the diblock copolymer systems. The 
dominant defects appear to be different between the sim- 
ulations and the experiment, leaving open the question 
of a general explanation for the coarsening behavior. 

For electroconvection (an anisotropic system), two 
main classes of patterns occur as the initial transition: 
normal and oblique rolls 0, 0, 0| . Electroconvection 
uses a nematic liquid crystal. The molecules of a liq- 
uid crystal align on average along a particular axis, re- 
ferred to as the director Normal rolls consist of a 
system of parallel rolls oriented with the wavevector par- 
allel to the director field (the average axis of alignment 
for the molecules). The main defect in the normal roll 
regime is isolated dislocations. The dynamics of these 
dislocations have been studied both experimentally and 
theoretically [H E EO, , though not in the context 
of domain growth. Oblique rolls correspond to stripes 
in which the wavevector forms a nonzero angle with re- 
spect to the undistorted director field. In the oblique roll 
regime, a quench typically produces a pattern consisting 
of domains of stripes with only two orientations, referred 
to as zig and zag rolls. In this case, the main defects 
are domain walls and dislocations. For a particular case 
of electroconvection, growth exponents of 1/4 were ob- 
served 0. However, in contrast to isotropic systems, 
disclinations were not present in this system. Again, this 
points to the need to better understand the dynamics of 
the various classes of defects if a general framework for 
understanding phase ordering in pattern forming systems 
is to emerge. 

In considering the existing work, there are at least two 
obvious questions. Given that the uniform and pattern 
forming systems exhibit similar topological defects, why 
are the resulting growth exponents so different? Given 
the range of growth exponents observed for pattern form- 
ing systems, can they even be explained in terms of the 
dynamics of topological defects? An important step in 
answering both of these question is the elucidation of the 
defect dynamics. The work with diblock copolymers has 
already made important contributions along these lines 
for disclinations. In this paper, we focus on the dynamics 
of confined dislocations found in the oblique roll regime 
of electroconvection plj . 

Figure^ is a top view of a pattern in an electrocon- 
vection cell that illustrates the defects of interest in this 
paper. The cell consists of a nematic liquid crystal con- 
fined between specially treated glass plates that align the 
director parallel to the plates along a single axis. This 
axis is defined to be the x-axis, or horizontal direction. 
The y-axis, or vertical direction, is also parallel to the 
glass plates, but perpendicular to the undistorted direc- 
tor. The z-axis is taken perpendicular to the plates. The 
nematic liquid crystal is doped with ionic impurities. An 
ac voltage is applied perpendicular to the plates. There 
exists a critical value of the applied voltage, T4, at which 
a transition from a spatially uniform state to a striped 



state occurs. The striped state consists of convection 
rolls with a corresponding periodic variation of the direc- 
tor and charge density. As shown in Fig. ^ we studied 
the case of oblique rolls, and the two classes of topo- 
logical defects are domain walls and dislocations. The 
dislocations are special in that they are mostly confined 
to vertical domain walls. Disclinations do not occur, as 
the stripes are not easily curved. As we will show, the 
domain walls are essentially static, so only the dynamics 
of the dislocations are of interest. In this paper, we re- 
port on qualitative features of the dislocation dynamics 
and their interactions. In addition, we will report on the 
time dependance of three global measures of defect ki- 
netics: dislocation density, average domain wall length, 
and the total domain wall length per area. The rest of 
the paper is organized as follows. Section II describes 
the experimental details. Section III presents the results, 
and Sec. IV is a summary and discussion of the results. 



II. EXPERIMENTAL DETAILS 

The details of the experimental apparatus are de- 
scribed in Ref. The nematic liquid crystal N4 
was doped with 0.1 wt% of tetra n-butylammonium bro- 
mide [(C4Hg)4NBr]. Commercial cells j2j| with a quoted 
thickness of 23 /im and 1 cm^ electrodes were used, giv- 
ing an aspect ratio of 435. The average wavelength 
of the rolls was 51 /im. The sample temperature was 
maintained at 30.0 ± 0.002 °C. The patterns were ob- 
served from above using a modified shadowgraph setup 
[2^ 13 nil that emphasized the contrast between zig 
and zag rolls. The magnification was chosen to monitor 
the largest possible area of the sample while maintaining 
enough resolution to resolve the stripe pattern. The area 
imaged contained approximately 150 rolls. 

All of the results reported here are for a fixed quench 
depth of e = /V^ — 1 = 0.05. After a quench, domains 
of zig and zag rolls form within 60 seconds. (The rele- 
vant relaxation times for electroconvection are the direc- 
tor relaxation time (0.6 s) and the charge relaxation time 
(4 X 10~^ s).) The domains are separated by two classes 
of grain boundaries. There are domain walls that extend 
horizontally. Across these grain boundaries, the phase of 
the stripes is continuous. The vertical grain boundaries 
consist of discrete steps that are formed by dislocations. 
Each step contains either an extra zig or zag roll. Both 
of the grain boundaries are illustrated in Fig.^ schemat- 
ically in Fig. ^ and using an actual image of the system 
in Fig. Figure^ illustrates the Burgers vector con- 
struction for one of the dislocations, demonstrating the 
extra 27r of phase that occurs when traversing a closed 
path around the dislocation core. Vertical boundaries on 
opposite sides of a domain are always composed of dislo- 
cations of opposite sign. Also shown is the definition of 
the domain wall length, Li, used later in the analysis of 
coarsening. 

The dislocation number density, p{t) = n{t)/A, was 
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measured by counting the number of dislocations {n{t)) 
in a fixed viewing area A. Because of the resolution used 
to take the images, and the fact that the defects were 
in such close proximity to each other within the domain 
walls, the defects were counted by hand. At each time, 
p{t) was averaged over ten different quenches. For the 
same times, we measured the length of the individual 
horizontal domain walls in the system, Li (see Fig. ^). 
From this, we computed both the average horizontal wall 
length (< L >— jj'^Li, where N is the number of 
horizontal walls) and the total domain wall length per 
area {L = Both of these quantities were av- 

eraged over seven quenches. The results are plotted in 
Fig. IHl Here time is scaled by the director relaxation 
time, Td — 0.6 s T^. In order to avoid issues of the ini- 
tial growth of the amplitude and wavenumber, we only 
considered times greater than AOOrd- After this time, the 
average wavenumber changed by less than 0.3%. 



III. RESULTS 

The first obvious difference between this system and 
previous experiments is the dramatic anisotropy of 

the growth. The dislocations in the vertical grain bound- 
aries move essentially horizontally. Some discrete steps 
in the vertical direction are observed. This is very differ- 
ent from the usual glide or climb of a dislocations. This 
is best seen in a movie of the motion. An archived movie 
made from snapshots taken every 60 seconds and with a 
playback rate of 0.3 s/frame can be found at Ref. [26l |. 
This horizontal motion is directly responsible for changes 
in the horizontal length scale of the domains. 

In contrast to the dislocations, the horizontal grain 
boundaries are effectively stationary. Any vertical "mo- 
tion" of domain boundaries occurs when an individual 
horizontal grain boundary is eliminated. This elimina- 
tion occurs when the oppositely charged dislocations that 
form the domain wall's endpoints(see Fig. annihilate. 
If only a single pair of endpoints annihilate, one has a dis- 
crete change of the vertical length by 49 microns. This 
step size is set by the step height corresponding to one 
dislocation in a vertically orientated domain wall. In or- 
der to eliminate a domain completely, the entire set of op- 
positely changed dislocations that form the vertical walls 
must annihilate each other. This is illustrated by the se- 
ries of snapshots in Fig.|21taken from the archived movie 
[2^ . In this case, the number of dislocations on each 
side of the domain differ by one. After the elimination of 
the domain, a single dislocation remains. Many domains 
contained equal numbers of dislocations, producing no 
isolated dislocations. 

There are a small number of isolated dislocations, as 
seen in Fig. |21 These dislocations move through the sam- 
ple either by climbing or gliding. Isolated dislocation ei- 
ther eventually move into a domain wall or annihilate 
with an oppositely charged dislocation in a domain wall. 
Figure, iniis a series of snapshots taken from the archived 



movie '2q that illustrate the incorporation of an isolated 
dislocation into a domain wall. This is consistent with 
the behavior of isolated dislocation predicted in Ref. [U . 
They rarely annihilate with other isolated dislocation be- 
cause of their extremely low density. 

During the evolution process, the dislocations often 
exhibit behavior substantially more complex than that 
shown in Fig.|2| For example, Fig.^Ji-d illustrates a case 
in which a single dislocation in a domain wall is observed 
to be undergoing oscillatory motion. Eventually, this de- 
fect is close enough to an oppositely charged defect that 
it accelerates toward that defect and is annihilated (Fig. 
0^ - h). Figure illustrates another interesting behavior. 
(Th is event is can also be found in the archived movie 
|2g.) Entire walls of oppositely charged defects can ap- 
proach within some distance, and then move apart. In 
a simple picture, oppositely charged defects would at- 
tract each other and such a bounce would not be possible. 
However, these defects are not isolated, and one must ac- 
count for the entire field of defects to describe these more 
complicated motions. This is outside the scope of this 
paper and will be the subject of future work. One also 
observes initially coherent walls of dislocations undergo- 
ing dispersion as they move (Fig.EI). This separation of 
dislocations occurs because the details of the interactions 
result in dislocations within a wall moving with different 
speeds. Finally, an extremely rare event is the nucleation 
of a new zig or zag domain. What is interesting is that 
extremely thin domains are possible in which the hori- 
zontal walls are slightly curved and the ends to contain 
dislocations (see Fig. Ejl. 

Given the complex nature of the dislocation dynamics, 
for the purposes of this paper, quantitative measures fo- 
cused on the time dependence of global quantities. The 
three quantities of interest are the dislocation density, the 
domain wall density, and the total domain wall length. 
The time dependence of these quantities is able to pro- 
vide information about the rate of ordering along the 
axes parallel and perpendicular to the director. The re- 
sults are summarized in Fig. |S1 

One observes that all of the quantities are consistent 
with power-law growth after t ~ lOOOr^. Fits of the data 
in this regime give: p ~ t-0-32±o.02. ^ _ ^-o.i8±o.o7. 
and < L >~ iO-33±o.02 -^j^jj ^j^g limited range that we 
are able to observe, it is difficult to show conclusively 
that a scaling regime has been reached. However, two 
conclusions are clear: (a) the growth is anisotropic and 
is best described by at least two different exponents; and 
(b) the exponents are different from previously studied 
striped systems 0, 0| . The first point is illustrated in 
Fig. IHfc, where the dashed line has a slope of —1/3. The 
second point is illustrated by the dashed curves in Fig.|S^ 
and b. Here, the dashed curves represent the previously 
observed growth exponent of 1/4. 

Given that the defect motion is essentially confined 
to the horizontal (or x) direction, we tested for the ex- 
istence of two different growth exponents. Typically, 
for domain growth, the scaling hypothesis assumes all 
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lengths are scaled by a single scale factor R{t) ~ t", 
and area scales as R^. For our system, because of the 
obvious anisotropy, we assume that horizontal lengths 
and vertical lengths of domains are scaled independently 
by scale factors ^ t"^ and Ly ~ respectively. 
Therefore, the number of domains, c, in our fixed view- 
ing area A scales as c ^ A/{LxLy). Because the dis- 
locations are confined in vertical walls, the number of 
defects is given by the number of domains times a typi- 
cal vertical dimension of the domain: n{t) ~ cLy. This 
gives pit) = nit) /A - Ly/{L^Ly) ^ t""-. The aver- 
age horizontal domain wall length (< L >) scales as 
< L >~ Lx ^ As with the number of defects, 

the total horizontal domain wall length scales as cLx- 
Therefore, the total domain wall length per viewing area 
will scale as L ~ cLx/A = L^/iL^Ly) ^ The re- 

sults in Fig.|Slare consistent with two growth exponents, 
Ux — 1/3 and ny — 1/5, with the growth in the vertical 
direction being substantially slower than the growth in 
the horizontal direction. Here, the time dependence of 
p{t) and < L > provide independent measures of Ux- 

Another issue is the impact of the finite viewing win- 
dow, which determines the maximum time for which we 
can view the system. This is seen in Fig. (SJ;, where the 
possible scaling of L{t) breaks down for t > SOOOr^, even 
though the other measures still exhibit possible scaling. 
This is understandable because this is the time when 
most horizontal boundaries extend across our field of 
view, and yet contain significant numbers of step dislo- 
cations of opposite sign. Therefore, as these dislocations 
annihilate with each other, both p(t) and < L > continue 
to evolve, while L(t) remains effectively constant. 

IV. SUMMARY 

We report observations of the dynamics of an inter- 
esting class of dislocations: dislocations confined to ver- 
tically oriented domain walls between two degenerate 
oblique rolls. The dislocations exhibit highly non-trivial 
behavior, forming coherent grain boundaries, exhibiting 
motion that is neither climb nor glide, and executing in- 
teresting dynamics, including oscillations. The motion of 
isolated dislocations toward the horizontal grain bound- 
aries and the confinement of the dislocations to grain 
boundaries is expected on general arguments from ampli- 
tude equations |2l| . However, existing theoretical work 
focused on the horizontal domain walls and individual 
dislocations |2^. A more detailed theoretical and exper- 
imental study is required to understand the full range 
of observed behavior. For example, it is clear that the 
observed dynamics are often due to the interactions of 



many dislocations. 

Measures of the global properties of the topological de- 
fects provide some insight into the phase ordering of this 
system. On the one hand, the combination of anisotropy 
and stripes results in a relatively simple system. There 
are only two classes of defects (domain walls and dis- 
locations), and their basic motions are straightforward. 
The dislocations move horizontally, and the domain walls 
are essentially stationary. This suggests that the phase 
ordering should be relatively straightforward to under- 
stand. On the other hand, the system is an interesting 
example of how stripes can make the system more com- 
plicated than the standard uniform systems. This is best 
seen by comparing the system to two standard universal- 
ity classes for phase ordering in uniform systems: Ising 
and X-Y models. 

Aspects of the electroconvection patterns are analo- 
gous to an Ising system, i.e. a systems of spins with two 
states. In our system, the two states are the zig and zag 
rolls, and the horizontal domain walls are the topological 
defects one would expect in an Ising system. However, 
the presence of the stripes results in additional topolog- 
ical defects: dislocations. These defects are analogous 
to vortices in an x-y model, i.e. spins with any orienta- 
tion in the plane, in that both vortices and dislocations 
have the same topological charge. In our system, they 
occur predominantly in vertical domain walls or as steps 
between two horizontal domain walls. Even with this 
ambiguity, because there is no obvious conservation law 
(both zig rolls and zag rolls are eliminated), one expects a 
growth exponent of 1/2 for both the Ising and X-Y model 
12 . The observed growth for electroconvection is clearly 
slower, being consistent with 1/3 in one direction and 1/5 
in the other. This suggests that the confined dislocations 
represent a new type of coarsening dynamics. 

One direction that may provide further insight into the 
exact cause of the slower dynamics is comparisons with 
simulations of an anisotropic Swift-Hohenberg model 
|27j |. Initial simulations are consistent with a growth ex- 
ponent of 1/3 for shallow quenches with pinning effects 
becoming important as a function of quench depth 27]. 
This will be the subject of future work. 
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FIG. 1: (a) A schematic drawing illustrating a transition from 
a zig to a zag and back to a zig region of tlie system. The black 
lines represent lines of constant phase between the stripes. 
There are two horizontal domain walls, one at each transi- 
tion. Also, one example of each type of step dislocation is 
shown. The Burgers vector is also illustrated for one of the 
dislocations. The length of a horizontal domain wall is in- 
dicated by Li. (b) A close up of a section of the system 
that shows two step dislocations (black squares). The white 
squares are provided as an aid in counting the zig and zag 
rolls. One finds that there are two more zig rolls than zag, 
as expected with two dislocations. The dashed line highlights 
one of the horizontal domain walls. The scale bar represents 
0.15 mm. 

FIG. 2: Four images separated by 60 seconds illustrating the 
collapse of a zig domain. The scale bar represents 0.15 mm. 

FIG. 3: Images separated by 150 seconds illustrating the 
incorporation of an isolated dislocation into a domain wall. 
The scale bar represents 0.15 mm. 

FIG. 4: Eight images separated by 300 seconds illustrating the 
motion of a single dislocation highlighted by a circle in (a). 
A dust particle on the surface of the sample is not removed 
to provide a frame of reference. In (a) - (d), the dislocation 
oscillates. In (e) - (h), the dislocation is attracted to one of 
opposite sign and annihilated. The scale bar represents 0.15 
mm. 

FIG. 5: Images separated by 150 seconds illustrating oppo- 
sitely charge dislocations moving together, "bouncing", and 
ultimately attracting each other. One can also observe the 
spread of dislocations in a vertical wall that is highlighted in 
Fig.|S| The scale bar represents 0.15 mm. 

FIG. 6: Images separated by 300 seconds illustrating dislo- 
cations in a domain wall spreading out in time. The scale bar 
represents 0.15 mm. 

FIG. 7: Images separated by 150 seconds illustrating the 
nucleation of narrow domain without a dislocation pair. The 
scale bar represents 0.15 mm. 



FIG. 8: (a) Plot of log]^Q(p) versus log]^Q(t). Symbols are 
experimental data. The solid line has a slope —0.32. For 
comparison, the dashed line has a slope of —0.25. (b) Plot of 
logiQ(< L >) versus logio(^)- Symbols are experimental data. 
The solid line has a slope 0.33. For comparison, the dashed 
line has a slope of 0.25. (c) Plot of logjQ(I/) versus logjQ(t). 
Symbols are experimental data. The solid straight line has 
slope of -0.2. For comparison, the dashed line in (c) has a 
slope of -0.33. 
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